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We derive the continuum and lattice tree-level moving NRQCD (mNRQCD) through order mNRQCD 
is a generalization of NRQCD for dealing with hadrons with nonzero velocity The quark's total momentum 
is written as P'' = Mu*^ -\- k'^ where k'^ <^ Mu'^ is discretized and Mm^ is treated exactly. Radiative corrections 
to couplings on the lattice are discussed. mNRQCD is particularly useful for calculating B ^ tt and B D form 
factors since errors are similar at low and high recoil. 



1. INTRODUCTION 

Decays of heavy quarks are essential for deter- 
mining the CKM elements. Because of the mass 
difference between the initial and final states, 
these decays can lead to large recoil momenta. 
Current lattice simulations require that the in- 
verse lattice spacing, a~^, be greater then the 
largest momenta to keep errors reasonable. For 
the important case B ■nlvi, needed for V^b, the 
high recoil limit gives 2.5 GeV. This implies 

2> 2.5 GeV, which is much more then current 
simulations can handle. 

A solution to this problem has two parts. First, 
we use a moving B meson and so share the mo- 
mentum between the B and the tt. This allows us 
to find a frame in which all light quark momenta 
are the same order and errors are minimized. Sec- 
ond, defining the momentum of the heavy quark 
as 



where is the 4- velocity of the B meson, means 
that A:^ ~ (;^) Mu'^ < Mu^^. We then discretize 
/c^ while treating M^u^ exactly. This removes the 
large, uninteresting part of the 6's momentum and 
allows larger lattice spacings to be used. This pro- 
cedure is analogous to NRQCD where the heavy, 
stationary quark has momentum = M + fc^ 
and is discretized. 



2. DECAY KINEMATICS 

The B meson momentum, Pb ~ {P, / P) , 
can be split between the b quark and the specta- 
tor quark in its wavefunction: 




P„ - xPe 



Pb = {1-x)Pb -p±, 



where the choice of P describes the choice of 
frame. 

The light quark has transverse momentum 
p± ^ A ^ 500 MeV with Ap± ~ A, and its mo- 
mentum fraction is peaked at x 

Ax ^ Consequently, 



-jg- with range 



Au 



(2) 



which implies that discretization errors associated 
with the heavy quark match those of the light 
quark. Note also that ub ■ k ~ A even as u^, 
fc^ oo. 

Our goal is to find a frame, for the simulation, 
where the discretization errors are the same for all 
quarks. To see how the choice of frame affects our 
calculation we give three examples for B — > irlu 
in the high-recoil limit. In the P's rest frame, 
P = Mb, and 



P. 



k^Pu^K => P. >P„. (3) 



With this choice discretization errors associated 
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with the pion are far larger than other discretiza- 
tion errors. In the tt rest frame, P — mjr, and 



AMb 
2A'U 



P^ « Pu] (4) 



the i?'s discretization errors dominate. The best 
frame is the "spectator- tt Breit" frame, P = 
\/AMb, where 



TT c\ c\ 



(5) 



^ - P„ - fc - 800 MeV 

In this frame the errors are distributed evenly be- 
tween both hadrons, and an inverse lattice spac- 
ing ^ 2 GeV is adequate. This is a savings of 
j^j^ ~ i compared to the B rest frame making 
simulations 1000 times faster. 



3. mNRQCD DERIVATION 

To derive the tree-level effective Lagrangian for 
a moving b quark we use the FWT method to 
transform the QCD Lagrangian order by order in 
1/m into non- interacting quark and anti-quark 
pieces. Working through order 1/to^, this gives 
the well known HQET result 



C = 



Im 'mi 



a 



_ig_ 



A3 



Q (6) 



and the original wave function ij: can be written 
in terms of Q. 

To implement this formalism efficiently on a 
Euclidean lattice, we must convert it to an initial 
value problem (in time). To do this we first re- 
duce the above equations to a 2-component form 
and then remove higher order time derivatives us- 
ing a field redefinition. 

To recast the Lagrangian into 2-component 
form we write Q as a rotation of a spinor with 



two zero components: 



Q = A 



1 + 7 ■v)\( X 



^2 (1 + 7) \l{^-v) 1 + 7 







(7) 



Higher order time derivatives are removed using 
a field redefinition Afr, 



The resulting mNRQCD to order {K/mnf is: 

£ = (^Dt + iv- D 
1 



(8) 



(9) 



+ 



2m7 
1 

i 

(2 - 



- v-D 



D-E 



V 



rV ■ 



a- B 

Dx B 



+ 



i 



2 ^ {v-D,&^-2(v-D 
V- D.v- E 



a- Dx E' - E' X D 



Srn^ (7 - 



—a ■ i^v-D,vx E'^ 



O 



where m'' = (7, jv) and E' and B' are the electric 
and magnetic fields in the b quark's rest frame: 



E' = -f[E + vxB- 



7 



7 + 



-v(v.E)y (10) 



B'^j(B-vxE 



7+1 



V \ V ■ B 



(11) 



mNRQCD reduces to NRQCD when v ^ 0. 
Lower order versions of this formalism are con- 
sidered in |ll|,|. 

Currents and other operators involving the b 
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quark can be designed using the b quark field: 



5. Conclusion 



^27(1+7) 



(12) 



X 1 + 



2m 



8m' 



{2u- D I/>^- I/>1^ 



xAfr A 



O 



X 



Note that our asymmetric treatment of time 
and space means that our Lagrangian relies upon 
an expansion in m = rather then A/m as in 
HQET. Since \u\ ^ 00 as w ^ 1 this could be 
problematic. Our choice of frame for B —^ irlv, 
however, means that "fv in never larger then 0.35 
and the truncation errors are no larger then the 
l/nv' terms neglected in (^). 

4. Lattice mNRQCD 

We reformulate the mNRQCD Lagrangian 
of (H) on the lattice using the usual forms for 
lattice derivatives and fields as in j|] . To achieve 
high precision, we must include perturbative cor- 
rections to our tree-level coefficients. Renormal- 
izations through order are needed for the lead- 
ing order v ■ D: 



The method shown here is easily extensible to 
other processes with large recoil including B 
Dlv which is used to find Vcb, D — s- -kIv used for 
Vcd and B K*j Using mNRQCD, errors are 
comparable at at both low and high recoil. We 
can vary the momentum transfer continuously by 
varying the B meson velocity, ub- 
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X 



a 



(13) 



where ^ = 1 -I- ^las + £,2a1 

For D mesons, where -p- ~ i, we need to in- 
clude terms of order through tree-level and 
^ through order ag. This requires 5 additional 
renomalizations 

X^ [D^ - (v-dY +ja-B' + C3-fa-B] x 



X^ [ciD'-C2(v-D 



+X — C470' ■ V X E + C5- 



r 



■{,.b) 



(14) 



X 



7+1 

where Ci = 1 + Ci^iag. For B mesons, where 
■p- ^ jlj, we need ^ only through tree-level. The 
calculation of the necessary perturbative coeffi- 
cients is currently in progress. 



